The two-level atom laser: analytical results and the laser transition. 
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The problem of the two-level atom laser is studied analytically. The steady-state solution is 
expressed as a continued fraction, and allows for accurate approximation by rational functions. 
Moreover, we show that the abrupt change observed in the pump dependence of the steady-state 
population is directly connected with the transition to the lasing regime. The condition for a sharp 
transition to Poissonian statistics is expressed as a scaling limit of vanishing cavity loss and light- 
matter coupling, K ^ 0, <; — >■ 0, such that g^/K stays finite and g^/K > 27, where 7 is the rate of 
atomic losses. The same scaling procedure is also shown to describe a similar change to Poisson 
distribution in the Scully-Lamb laser model too, suggesting that the 1ow-k, \ow-g asymptotics is of 
a more general significance for the laser transition. 
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I. INTRODUCTION 

A single two-level emitter in interaction with a cavity 
mode is the simplest illustration of cavity quantum elec- 
trodynamics. Actually, in experimental situations the 
emitter, either a true atom or an artificial one (quantum 
dot), usually brings into the picture more than just two 
levels. Besides the pair of states which couple with the 
laser mode, other states are involved in the carrier ex- 
citation and scattering processes. It was early realized 
[l[ that incoherent pumping directly from the lasing lev- 
els destroys the coherence between these and the photon 
output is strongly reduced for large pump rates. This 
quenching effect acts against lasing and can be avoided 
by involving additional states. More recent realizations 
of single-emitter lasers 0, Q were modelled using at least 
four levels In a semiconductor quantum dot (QD) 

one typically encounters many-level, many-carrier situa- 
tions which give rise to multiexcitonic effects. Moreover, 
in the presence of a spectral continuum (e.g. wetting- 
layer extended states) where, in typical experimental sit- 
uations carriers are initially created by optical or electri- 
cal pumping, the creation of fermionic reservoirs becomes 
possible, from which carriers are injected into the QD 7]. 
The same charge distributions in the continuum are also 
a source of screening and additional dephasing @. All 
these effects are clearly beyond the reach of a two-level 
description. 

Despite this, the two-level emitter model still enjoys a 
wide popularity due to its simplicity and transparence 
and remains a subject of intense investigations Q. 
Also, models with more levels are sometimes mapped 
into effective two-level ones [l^, Ell- Many aspects of 
the light-matter interaction, like spectral properties as- 
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sociated with strong cou plin g, are quite well described in 
this framework 

An appealing feature of the two-level model is the fact 
that the carrier degrees of freedom are easily eliminated 
leading to closed equations for the purely photonic 
statistics. This is a feature shared with the Scully-Lamb 
random injection model , which addresses directly the 
photonic density matrix. Therefore these cases are par- 
ticularly appropriate for the study of the laser transition, 
and indeed, the Scully-Lamb model became a paradigm 
and a textbook case in this respect [l8l - l2l| . In compari- 
son, with a few exceptions [22, 23|, the two-level emitter 
has not enjoyed the same attention from this point of 
view. 

In the Scully-Lamb model it is shown that the tran- 
sition to lasing corresponds to the peak of the photon 
number distribution getting detached from the origin 
and moving away to larger and larger values. A ther- 
mal, exponential distribution changes into a Poissonian 
statistics. In terms of the Glauber-Sudarshan representa- 
tion [2^ , this means that the so-called quasi-distribution 
function, initially located around the origin of the com- 
plex plane, becomes narrowly concentrated on a ring of 
increasing radius. Such a behaviour was also obtained 
numerically [l|, |25'] or by using analytic approximations 
[2^ . in several different contexts. 

In the case of the two-level model, the numerical 
simulations show an abrupt change of regime, most 
clearly seen on the steady-state population behaviour: 
the strictly concave dependence of the upper-level pop- 
ulation on the pumping rate becomes with a high ac- 
curacy linear, in the interval of intermediate pump rates 
[lOl [Ta . [23 | . Various numerical evidence suggest that this 
transition is associated with the onset of lasing, but the 
connection to the behaviour of the Glauber-Sudarshan 
function, like in the random injection model, is far from 
obvious. 
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In the present paper we address the two-level laser 
problem by analytic methods. We show that the steady- 
state solution, usually obtained by the long-time limit of 
the Liouville-von Neumann evolution, can be expressed 
directly as a continued fraction involving the model pa- 
rameters. Exact relations as well as simple, yet accurate 
approximations for the solution can be derived in some 
cases from the full formula. The theory of continued frac- 
tions also provides error control for these approximations. 

But what is more important, the method also allows 
for a detailed characterization of the laser transition. The 
appearance of a jump in the solution (or some derivative) 
has its analog in the problem of phase transitions, to 
which lasing was sometimes compared, and the role of 
exactly soluble models is to clarify the mechanism of the 
jump. 

The two-level atom turns out to be such an exactly 
soluble, yet nontrivial example on which the laser tran- 
sition can be studied. The condition for the appearance 
of the above-mentioned abrupt change in the pump de- 
pendence of the population is proven to be equivalent to 
the vanishing of the Glauber-Sudarshan function at the 
origin. Moreover, one can show that this vanishing re- 
quirement is met only in a certain scaling limit for the 
model parameters, and then it also entails a Poissonian 
photon statistics. 

Similarities with the thermodynamic limit in phase 
transitions and appropriate scaling asymptotics, like the 
Grad limit or the weak-coupling, long-time limit, in the 
theory of markovian kinetics [26| come to mind. In all 
these examples a regime change occurs, rigorously speak- 
ing, only in a certain asymptotic parameter domain, and 
this is precisely what takes place in our case too. Of 
course, with parameters close but slightly away from the 
scaling limit the transition may still be seen, albeit less 
sharply. 

The paper is organized as follows: After defining the 
model (Secini the three-term recursion equation obeyed 
by an infinite sequence of photon expectation values [l6[ 
is discussed and its solution expressed as a continued frac- 
tion. It is shown that outside the interval of linear be- 
haviour the continued fraction converges rapidly and its 
truncation provides simple and accurate approximants 
(Sec. IIII[) . The discussion of the linear regime is more 
complicated and requires an extension of the recursion 
relation. This is achieved by noting that the photonic ex- 
pectation values are moments of the Glauber-Sudarshan 
function, and it is shown that the linear behaviour is 
equivalent to the condition that this function vanishes at 
the origin fSec. IIV|) . 

The next sections deal with the circumstances in which 
this condition is met. Heuristic suggestions in this 
sense are obtained using the rate-equation approxima- 
tion (Sec. El . Turning back to the full problem (Sec. IVl)) 
it is shown that exact vanishing occurs as a limiting case, 
in which the rate of cavity losses n and the light-matter 
coupling strength g go simultaneously to zero so that 
g'^ /k = const. A similar scaling limit was obtained by 



Rice and Carmichael, as a condition for a sharp lasing 
transition in the rate-equation approach to the many- 
emitter laser It should be noted that in contradis- 
tinction to this many-emitter result, in the present case, 
as the pump increases the system undergoes two transi- 
tions, a first one at the onset of lasing and a second one 
at larger pump values, when the lasing is destroyed by 
quenching. It is also shown that in this limiting proce- 
dure, in the lasing regime the Glauber-Sudarshan func- 
tion not only vanishes at the origin but also becomes 
(5-like distributed on a ring centered at the origin, in 
agreement with the expected Poissonian photon statis- 
tics. The radius of this ring starts increasing from zero 
when the pump reaches the threshold value but eventu- 
ally decreases back to zero due to quenching. 

The conclusions regarding the scaling limit and the 
laser transition go beyond the particular case of the two- 
level laser, as checked on the Scully-Lamb model (see 
Appendix |B]) . 

II. THE MODEL AND ITS STEADY-STATE 
BEHAVIOUR 

We consider a two-level emitter (atom or QD) in 
Jaynes-Cummings (JC) interaction with a cavity mode. 
Cavity losses, spontaneous recombination into non-lasing 
modes and incoherent pumping are included in the den- 
sity operator evolution through Lindblad terms with 
rates k, 7 and P, respectively. Assuming perfect reso- 
nance between the emitter and the cavity mode one has, 
in the interaction representation, the following Liouville- 
von Neumann equation (in h — 1 units) 

-p ^ ^z[Hjc, p] + ^ [2b pb^ - p - pb^b] 

+ [2v^c pc'v — c'cp — p c^c] 
p 

+ — [2c^v pv'^c — v^v p — pv^v'\ (1) 

The corresponding interaction Hamiltonian is given by 

Hjc = g [b^'v^'c + bcK] . (2) 

The creation and annihilation operators arc b^ , b for the 
laser mode, c\c for the upper (conduction band) level 
and for the lower (valence band) one. Only one 

carrier is present in the system, so that one has c'c + 
v'^v — 1. We use QD notations throughout, the atomic 
pseudospin formalism being easily recovered via a = v^c, 
(T^ ~ c^v, (T^cr = c^c, a.s.o. 

The Liouville-von Neumann equation can be solved nu- 
merically until the steady-state solution is reached. Then 
various steady-state expectation values, regarding level 
populations and photon statistics can be obtained as a 
function of the pump rate P or other parameters 

Examining families of such plots [8[ one clearly de- 
tects two types of behaviour. One example of each of 
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FIG. 1. Steady state solutions of the Liouville-von Neumann equation vs. pump rate P. Tiie parameters are: Set A (red lines): 
K = 0.01/ps, g — 0.1/ps, 7 — 0.02/ps and Set B (green lines) : k — 0.02/ps, g = 0.01/ps, 7 = 0.01/ps. (a) the upper state 
population Tie, semi-log plot. Circles correspond to the solution obtained by continued fraction (see Sec IIII[) . (b) the same in 
linear plot, showing the appearance of the linear middle regime, (c) average photon number n. (d) second-order autocorrelation 
function. 



these types is shown in Fig.l (see the Caption for the 
parameters used). The two cases are best distinguished 
by their pump dependence of the upper-level population 
Uc = {c'c). In a semi-logarithmic plot one may have, 
depending on the parameters, either a simple S-shaped 
curve or a double S-shaped one (Fig. la). The difference 
is more striking in a linear plot (Fig. lb), where in the 
first case one has a strictly concave function of P, while 
in the second there appears an almost perfectly linear 
shortcut in the middle P-region. This is strongly remi- 
niscent of the behaviour of thermodynamic potentials in 
the presence of a phase transition ^] . 

In the present case too, the transition from concave 
to linear dependence in the middle P interval is a clear 
simptom of an abrupt regime change in the functioning 
of the two-level laser. The figure shows that this "lin- 
ear middle" regime (as it will be subsequently called) is 
associated with inversion of population (ric > 1/2), a 
large photon output n — (Fig.lc), and a value of 

the second-order photon autocorrelation function at zero 
delay g(2)(o) = b'^ bb) / {h'^hf close to unity (Fig.ld). 
All these are indicative of the onset of lasing which, de- 
pending on the parameters, may take place or not, as 
illustrated by the two types of behaviour of Fig.l. 

In what follows it will be proven that the linear pump 



dependence of the population is indeed associated with a 
Poissonian photon statistics, characteristic for coherent 
light. At the same time, the conditions for the presence 
or absence of the transition, as well as cases in which the 
transition becomes sharp will be identified. 

III. SOLUTION BY CONTINUED FRACTIONS 

It is easy to see that the Liouville-von Neumann equa- 
tion leads to a closed system of equations for a subset 
of density matrix elements. In the basis of QD-states 
c) , \ v) and photon number n, these are the diagonal ones 
P'n,m Pn',n ^"^^ ^"^Y thosc off-diagoual elements between 
states having the same number of excitations c^c -I- b%, 
i.e. p™n+i their complex conjugates p^+i „. The lat- 
ter restriction is a consequence of the JC Hamiltonian, 
which conserves the excitation number. 

Equivalently, one may limit oneself to the following set 
of expectation values 

Cfc = (etc 6^=6^) , 

Vk = {v^vh'^H^) , fc = 0,l,2... and 

tk = -ig{v'^ch'^%^^^) , fc=l,2,3... (3) 
Obviously, one has Uc — (c^c) = cq and Uy — (v'^v) = wq. 
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and Co + t^o = 1. Of particular interest are the purely 
photonic expectation values 



Pk 



(4) 



and especially the low-index ones pi = n, and p2 related 
to the autocorrelation function by p2 = s'^HO) Pi- 

With the imaginary unit as a pref actor, the multi- 
photon-assisted polarization tk is real. From Eq. ([l} the 
equations of motion for these expectation values can be 
derived, together with the steady-state conditions 



d_ 

dt 
d_ 

dt 
d_ 

dt 



Ck 



Vk 



-{kn + "f)ck + Pvk - 2tk+i = , 



(5a) 



7Cfc - (kK + P)vk + 2ifc+i + 2ktk = , (5b) 



tk = g^ck + g^kck~i - g'^vk 



{2k - 1)k + P + 7 



ife = 0. (5c) 



By adding Eqs. (f5a| and (j5b| one obtains the steady- 
state balance relation between the cavity losses and its 
feeding through the photon-assisted polarization 



Kpk^2tk , k^ 1 



(6) 



Using this condition, Eq. (|5ap for k — can be written 
as P Hy = "fUc + Kn, which leads to 



P — nn 
P + 7 



(7) 



This shows that the knowledge of the photon output n is 
sufficient for the determination of the level occupancies. 

More generally, from Eqs. (|5a)) and ()5bp for arbitrary 
fc, the unknowns Ck and Vk can be eliminated in favour 
of tk, and using again the balance condition Eq. ^ one 
obtains closed equations for the photonic expectation val- 
ues Pk, in the form of a three-term recursion relation, first 
obtained by Agarwal and Dutta Gupta 



is available, namely po = 1. On the other hand, it is 
known (29j that the three-term recursion problem can be 
addressed using continued fraction techniques [H, [30l |. 
To this end Eq. is rewritten as 



Tk + A/c 



rk-i 



, fc > 1 , 



(10) 



with rk = Pk+i/Pk , Afe = Bk/Ak and fik = Ck/Ak ■ 

As before, the initial condition tq = n is not known. 
On the other hand it turns out that the large k limit is 
easy to obtain. Indeed, according to their definitions all 
Tfc should be positive and from Eq. the positivity 
requirement for both Vk-i and leads to 



Tk-l ^ 



Bk 



(11) 



An examination of Eq. ^ shows that the upper bound 
goes to zero as fc — > cxd, and therefore limfe_j.oo fk = 0- 

With this 'initial' condition at fc = oo, Eq. (jlOp is 
solved by a backward propagation of the recursion re- 
lation, rk-i — ^J.k/i'^k + Tk)- This leads to a continued 
fraction expansion of the solution, in terms of the param- 
eters ^k a-iid Afc 



ro = 



Ml 



Ml 



Ai + ri fi2 
Ai + 



A2 + r2 
Ml 



(12) 



Ai + 



M2 



A2 + 



M3 



A3 



The continued fraction is convergent [3l| and the result 
is numerically in perfect agreement with the long time 
limit of the Liouville-von Neumann solution (see Fig. la). 

(k) 

Approximate algebraic expressions Tq for the solution 
ro are obtained by setting in Eq. = 0. Fig. 2a 

shows a comparison of the full solution n = with its 

(2) 



Ak Pk+i + BkPk — Ck Pk-i = , fc ^ 1 , (8) second approximant rg which is a rational function of P 

whose expression follows straightforwardly from Eq. Q 



with 



Ak = 



2k 



Bk = 



Ck 



kK + P + j ' 

fc K — P + 7 
+ 

kK + P + "f 



k K 



(fc - l)K- 
(2fc- 1) 



-p + 7 
K + P-t-7 



kP 

(fc - 1)k + P + 7 



(9) 



This is an important feature of the model, allowing direct 
access to the photon statistics after the elimination of the 
information regarding the atom subsystem. 

Solving Eq. ([5]) recursively requires knowledge of two 
initial conditions, and the difficult point is that only one 



and r^^^ = CiB2/{BiB2 + A1C2). The corresponding 
approximants for ric obtained by Eq. ([7]) are shown in 
linear (Fig. 2b) and semi-log (Fig. 2c) plots respectively. 
The agreement is remarkable on the whole range of pump 
values in the case when the transition is absent (Fig. 2d), 
and when the transition does take place, the agreement 
remains good for the pump values away from the linear 
middle interval. 

A closer inspection of the numerics shows that inside 
this interval the continued fraction becomes much more 
slowly convergent. The reason for this is the fact that 
the first values of Bk, and with them those of Xk, are 
negative, resulting in negative values for the low-ranking 
approximants of the positive quantity tq. Indeed, from 
Eq. ^ it is clear that Bi may take negative values due 
to its first term. For this to happen P should not be too 
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small. On the other hand for large values of P the last 
term becomes dominant, making Bi positive again. In 
other words, Bi can be negative but only for an appro- 
priate choice of the parameters k, 7 and 17 , and only for 
intermediate pump values. This makes the analysis of 
the linear middle interval look rather complicated. 

Surprisingly, the numerical evidence shows that the sit- 
uation is quite simple. One remarks that the full solution 
obeys to a good accuracy the approximation provided by 
the ansatz 



with ro as in Eq. (jl2[) and similarly ri given by 



''0 = -Ao 



Aq 2k 



8g^ 



(13) 

which leads, using Eq. ([7]), to a linear expression for the 
upper-state population 



1 k(P + 7-k) 



(14) 



Moreover, in terms of the recursion relation, this ansatz 
is equivalent to the assumption that Eq. would hold 
even for k = with the last term CoP-i supposedly 
vanishing. Such an interpretation was put forward by del 
Valle and Laussy [HI, based on the argument that Co = 
0. But this definitely requires a more careful analysis, 
because on the one hand Co comes multiplied by the ill- 
defined p-i and on the other hand the ansatz is obviously 
not working for all parameter values. 

Nevertheless, it is clear (see Fig. 2) that the values given 
by Eqs. and (fH)) describe successfully the linear 
middle regime. Moreover, since this regime is suppos- 
edly connected with lasing, it is important to clarify the 
deeper meaning of the ansatz. In other words one has to 
answer two basic questions: (i) if indeed Eq. ^ can be 
extended for fc = 0, which is the right interpretation of 
its last term in this case and (ii) why is this term van- 
ishing for parameters corresponding to the linear middle 
regime. 

The first question is answered in the next Section and 
the second in Section IVTl 

Before closing this Section two remarks are in order. 
First, it is easy to see that if Bi is positive, all Bk are pos- 
itive too, and then the continued fraction coefficients are 
positive. In this case the odd approximants Tq^'^^^' 
form a sequence of decreasing upper bounds for the true 
solution To, while the even approximants r^'^'' give rise to 
an increasing sequence of lower bounds. This provides an 
efficient control of the truncation error, when the system 
is away from the linear middle regime. 

Second, we show that the continued fraction solution 
provides exact analytic expressions for the autocorrela- 
tion function for both the low- and the large-pump limit. 
One has 



.(2)(o) = ^^!l 

Pi ro 



(15) 



ri 



M2 



A2 + 



A*3 



(16) 



A3 + 



As P all vanish linearly while stay finite. 
Therefore, in this limit one has 



and as a result one obtains 



(0)|p=o 



4g2 



K7 



K -I- 7 

3k + J Ag^ + k{k + 7) 



(17) 



(18) 



In the large P asymptotics /i/j ^ Pk/2K grow linearly 
with P and therefore are dominated by the quadratically 
increasing Afe ^ P^/8g^. In this case too the asymp- 
totic behaviour is given by Eq. (|17l) with the outcome 
hmp^oo 5''^'(0) — 2. These facts are borne out by the 
numerical results. 



IV. GLAUBER-SUDARSHAN 
REPRESENTATION 

An alternative approach to the photon statistics is 
based on the expansion of operators in the photonic 
Hilbert space as integrals over coherent-state projectors 
[2^ . This is known as the Glauber-Sudarshan (GS) or 
phase-space representation. 

In our case the Hilbert space is enlarged by the QD 
degrees of freedom and the density operator can be rep- 
resented as a 2 X 2 matrix made of blocks acting on the 
photonic degrees of freedom. For each of these blocks the 
GS expansion can be used 



Pc 



da 



Rc aQ 
a*Q* 



(19) 



The GS functions Rc and Ry representing p'^'^ and re- 
spectively p"", are in principle functions of the complex 
argument a. Nevertheless, as discussed in Sec. piip the 
elements of the density matrix which are diagonal in 
the QD-level index are diagonal in the photon number 
too, and this entails that the corresponding GS functions 
depend only on the radial parameter s — \a\^ and are 
angle- independent Rc = Rc{s), Rv — Rv{s). Similarly, 
the off-diagonal block p'^" has non-zero entries only one 
step above the diagonal. This corresponds to a GS func- 
tion with the structure aQ{s). Also, in accordance with 
Eq. ^ it is convenient to introduce T(s) — —igQ{s), 
which will turn out to be real. 

The advantage of using the GS representation becomes 
obvious after expressing the expectation values of Eq. ^ 
in this formalism. Indeed, the procedure is convenient 
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FIG. 2. Comparison between the full solution (red line), the second approximant to the continued fraction (green squares) and 
the linear middle ansatz, Eqs. (|13I14|I (blue circles). The parameters are the same as in Fig.l: Set A for panels (a), (b), (c), 
and Set B for panel (d). 



for calculating averages of normal-ordered photonic op- 
erators 20, 24]. One obtains 



Cfe = / Rc{s)s''ds , Vk= Rv{s)s''ds , 

/Q Jo 

f oo 



/•OO 

Jo 



(20) 



This shows that the quantities of interest are moments 
of the GS functions, and it is clear that this allows to 
extend their definition. For instance, now it is possible 
to define to according to Eq. (PU)) . whereas in Eq. ^ it 
made no sense. 

The purely photonic expectation values are obtained, 
as in Eq. (j4]), by tracing over the QD-indices 



POO 

Pk ^ Ck+Vk ^ / R{s)s''ds , 

^0 



(21) 



in which R{s) = Rc{s) + Rv{s). For further reference 
let us remind here that the photon number probability 
distribution pn n — P^n + Pn"n is given in terms of i?(s) 
by 



Pn,« = / R{s)—e "ds 
lo n\ 



(22) 



The Poissonian statistics is generated by a S-\ike GS func- 
tion, R{s) = 6{s — So), which corresponds in terms of 
the photonic expectation values pk to an exponential be- 
haviour Pk = So- 

The translation of the Liouville-von Neumann Eq. ([T]) 
into the phase-space language leads to Fokker-Planck 
equations for the GS functions [13], obtained by the cor- 
respondence 6 — > a and 6^ — >■ a* — d/da when the pho- 
tonic operators act on p from the left and b a — d/da* 
and 6^ — a* when they act from the right. Taking into 
account that the GS functions appearing here depend 
only on the radial argument, one obtains the counter- 
part of Eqs. dS]) in the form 



d_ 
dt 



Rr 



-2sT + k{Rc + sR'^) 

-jRc + PRv=0 , (23a) 
-2sT' + 2{s - 1)T + k{R., + si?;) 



+ -fRc-PRv=0, 
R,)~g^R', + K{^T + sT') 



T = 



(23b) 



(23c) 



where the prime stands for the derivative with respect to 



7 



s. By adding Eqs. (I23ap and (j23bp one obtains succe- 
sively 

- 2(sT)' + nisR)' = , skR = 2sT , (24) 

whcrefrom the analogue of the balance condition Eq. ([5]) 
follows 



kR = 2T . 



(25) 



The link to the equations of motion Eqs. ([S]) is ob- 
tained by taking the k-th moment of both sides of the 
Fokker-Planck Eqs. ([23]) and using integration by parts 
for the terms containing derivatives. The advantage of 
this approach is that not only the equations of motions 
are recovered, but now the k = version of the last one 
(Eq. ([5c|) ) can be obtained, and reads 



g^CQ-g^vo+g^RciO) 



-K + J + P 



, (26) 



which amounts to extending the interpretation of kck-i 
for fc = by the value of the GS function Rc{s) at the 
origin. The process of eliminating Ck and Vk in favour of 
tk can be repeated, but now including the k — Q case. 
The balance relation Eq. ([6|) holds for fc = too, owing 
to its GS version Eq. As a result one recovers the 

three-term recursion formula Eq. (|8]) extended for fc = 
as 



P 



-K + 'J + P 



R{0) = 



(27) 



This shows that RiO) is the correct interpretation of 
kpk-i for fc = [34I. In deriving the last equation we 
made use of Eq. (|23a[) which, for fc = and s — 0, 
provides the relationship between i?c(0) and i?^(0) — 
R{0) — RciO), allowing to express i?c(0) in terms of i?(0). 

Having established the generalization of the recursion 
relation to fc = 0, it becomes clear that the ansatz con- 
tained in Eq. ([T3| . and which leads to the linear middle 
behaviour Eq. is equivalent to the vanishing of i?(0). 
We will therefore be concerned in what follows with the 
circumstances in which this vanishing may take place. 

One way to address this problem is to eliminate Rc{s) 
and Ry{s) from Eqs. in order to set up a differen- 
tial equation for T(s) which, through Eq. (|25p . provides a 
closed equation for the purely photonic GS function R(s). 
This essentially parallels the procedure of obtaining the 
recursion relation, except that we perform it at the level 
of the GS functions and not of their moments. One ob- 
tains in this way a second-order differential equation of 
the form 



{as^ + bs^)R"{s) + {c + ds + es^)R{s) 

+ {u + vs + ws^)R{s) = 



(28) 



Details about the derivation of this equation and the ex- 
pression of the coefhcients in terms of the model parame- 
ters K, 7 and the pump rate P are given in the Appendix 

m 



The structure of this equation corresponds to s = 
being an irregular singular point [SSj]. Only one of the 
fundamental solutions behaves around s = as a power, 
R{s) ^ s'^ and the analysis leads to the result (7 = 0. 
The other fundamental solution is highly singular R{s) ^ 
exp{X/s), A > [l^l- These conclusions are parameter 
independent and show that one cannot find a solution 
which vanishes at the origin. 

This seems to contradict the previously mentioned nu- 
merical findings. The paradox comes from the require- 
ment i?(0) — 0, which is too strong. In fact the numerical 
evidence hints only at i?(0) being negligibly small. We 
show below that, indeed, the strong condition R{0) = 
holds only as a certain limiting case, which corresponds 
to a sharp laser transition. Away from this limit the 
transition gets smoother and the condition is expected 
to hold only as an approximation. This is a familiar 
situation, also encountered in connection with the many- 
emitter lasers The scaling limit procedure found in 
that context is based on an analysis performed at the 
level of the laser rate equation. In the present case too, 
we first take guidance from the rate equation approxi- 
mation to find the appropriate limiting regime, which we 
later analyze at the level of the exact solution. 



V. THE RATE EQUATION APPROXIMATION 

The infinite hierarchy of the equations of motion 
Eq. ^ can be interrupted by assuming that higher-order 
expectation values can be approximately factorized into 
products of lower-order ones. The factorization leads to 
the rate-equation approximation to the full Liouville-von 
Neumann problem |27| . In our case, using 



ci = (c^'cb'^b) « (c^c) (b^b) = copi , 



(29) 



and limiting ourselves to the steady state, we are left with 
the following closed set of equations for the quatities of 
interest ric = cq, n.v = vq and n = pi 



Uc + Uy = 1 , 
— 7 Tic + P Uy = Kn , 

T Uc {n + I) — T Uy n = K n . 



(30) 



Here F — Ag"^ /{n + 7 + P) is the rate of spontaneous 
emission into the laser mode, and the physical interpre- 
tation of different terms is obvious. Let us mention only 
the fact that the decrease of T with P is a reflex of the 
dephasing effect of the pump, which leads eventually to 
the quenching of photon output. 

Solving the first two equations for the level occupancies 
in terms of the photon number leads to 



P 



K n 



J + Kn 
1 + P ' 



(31) 
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which, introduced in the last one, gives rise to a quadratic 
equation for n 



2 (Kn)^+ 



4g2 
kP = , 



{k n) 
(32) 



The two solutions have opposite signs and, of course, only 
the positive one is physical. 

For the sake of the following discussion the above equa- 
tion was rewritten in terms of the rescaled photon num- 
ber n — Kn. In this form, in the limit At — > the free 
term disappears and one of the roots becomes n = 0. 
This vanishing solution holds in the parameter domain 
in which the other root is negative but, when positive, it 
is this second root that is retained. Therefore an abrupt 
change of behaviour takes place when the non-zero solu- 
tion changes sign. As a function of the pump rate the 
second root is a quadratic polynomial, initially increas- 
ing but then decreasing, as predicted by the quenching 
effect. This behaviour is preserved even in the k — )■ 
limit, provided one takes simultaneously 5 — in such a 
way that = j n stays finite. In this limit the solution 
reads 

" = ^ 8^ ^ "^^^ 

when ^{P) ^ and n — otherwise. For the upper-state 
population this leads to 



852 



if i/(P) < 
if i/(P) ^ 



(34) 



These solutions are shown in Fig. 3. It is clear that the 
concave dependence of ric on P is interrupted by a linear 
middle behaviour, in the interval defined by the roots of 

The situation described by the rate equation approxi- 
mation suggests that, in the above-defined scaling limit, 
the interval of positive values for i^(P) coincides with the 
interval in which i?(0) = 0, and in which the linear mid- 
dle ansatz Eq. (|13p holds. This is indeed the case, as will 
be proven in Sec. I VII 

Several comments are in order. First, the appearance 
of the transition is conditioned by existence of real, posi- 
tive roots of the equation J^(P) =0. It is easy to see that 
real roots exist if ^ 27, or g^ ^ 2/c7, and then they 
are both positive. The two types of behaviour discussed 
in Fig.l are distinguished by whether this condition is 
fulfilled or not. 

Second, the scaling procedure defined here k — 0, 
(7—^0, g^ / K finite, corresponds entirely to the limiting 



0, /?//t finite) considered by Rice and 
27| , because the /3-factor is proportional to 



regime (k — > 0, 
Carmichael 

g^. Moreover, the Jaynes-Cummings coupling constant 
appears as an explicit parameter in all laser models and 
therefore the scaling limit formulated in terms of g is not 
limited to the two- level atom model (see Appendix IB|) . 



Finally, since n is the rate of photon emission out of 
the cavity, its nonvanishing value as k goes to zero ex- 
presses the appearance of a finite laser output even as 
we get arbitrarily close to a perfect cavity, by a corre- 
sponding increase of the photon generation inside. This 
is physically interpreted as "an explosion of stimulated 
emission" [27|, typical for the lasing regime. 

The rate equation approach provides only partial in- 
formations about the photon system. The full photon 
statistics is contained in the GS function R{s), whose 
scaling behaviour is addressed in the next Section. 



VI. THE SCALING LIMIT OF R{s) 

As n — >■ c», the moments of R{s) cease to exist, show- 
ing that the GS function tends to move its weight on 
larger and larger s values. In this situation, meaningful 
quantities are the rescaled ones, like 

/•OO 

h = Pi = / KR[s)sds. (35) 
Jo 

In the rescaled variable t = ks, one has 

/•OO 

Pi = / R{t)tdt , (36) 



where the rescaled GS function R{t) is defined by 



Rit)^-R(-] . 



(37) 



Similarly, all the higher moments pk are rescaled as pk = 
K^Pk and become corresponding moments of R{t). 



Pk 



R{t) t'' dt 



(38) 



The task now is to establish the equation obeyed by 
R(t) in the scaling limit and to study the behaviour of 
its solution around the origin. Intuitively, according to 
Eq. (1571) . the graph of R{t) is obtained from that of R{s) 
by compressing the latter by a factor of 1/k along the 
abscissa and expanding it by the same factor along the 
ordinate. This would bring the rescaled function to S{t), 
in the limit k — s> 0, were it not for the opposite tendency 
of R{s) to move away from the origin, as discussed above. 
The net result of these competing trends is analyzed be- 
low. 

To this end, one translates Eq. ([25]) obeyed by R{s) 
into an equation for R{t), retaining in the resulting co- 
efficients of the latter only the dominant terms in the 
scaling parameter k. One obtains 34] 



^2 
4.g2 



,7 + ^^ , 
u{P)]R. 



I t 



1 



P 



kRI 



(39) 
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Pump (1/ps) Pump (1/ps) 



FIG. 3. (a) n (full red line) is equal to (blue dashed line) when the latter is positive and zero (green dashed line) otherwise, 
(b) the values of the upper level population Uc for these cases, Eq. ((34]), showing a linear middle behaviour. 



The appearance of the small parameter k with the deriva- 
tives suggests a WKB approach to the k ^ asymptotics 
of the solution. In other words one searches the solution, 
up to a normalization factor, in the form 



(40) 



in which Lp{t) is taken in the leading, zeroth order in k. It 
is clear that when k, gets smaller, the value of R{t) around 
the minimum of Lp[t) is greatly enhanced, at the expense 
of the values at other points which, in the view of nor- 
malization, become negligible. In the limit one obtains a 
(5- function concentrated at the minimum of ^{t). 

The equation obeyed by (p{t) in the leading order has 
the form of a quadratic equation for its derivative 



4 {^'f+ 



It 



1 



p 







(41) 

Around t — one of the roots behaves like if' ~ 2g^P/t^, 
i.e. if ~ —2g'^P/t which, by Eq. PO)) . leads to the sin- 
gular solution mentioned in Section IIVI The regular 
solution comes from the other root, for which ip'{0) — 
-2u[P)/P. 

Two cases arise, depending on the sign of v{P): (i) As 
long as v{P) is negative one has > and t ~ 

is a minimum for ip{t). Then, according to the above 
discussion, R{t) tends to 6{t) in the scaling limit, (ii) 
When z^(P) becomes positive, ip' starts aX t — with 
negative values and crosses the abscissa at t — ly. Then 
ip{t) has a local maximum at the origin and therefore 
the values of i?(0) become vanishingly small in the limit 
K — 7- 0. This is precisely the requirement for the linear 
middle behaviour to hold. The values of R{t) are now 
concentrated at the point of minimum t = v{P). In other 
words R{t) 5{t - v). 

In the first case, according to Eq. , all rescaled ex- 
pectation values are vanishing in the scaling limit, except 
^0 = 1; while in the second regime they become pk = ■ 
This exponential behaviour of the normal-ordered aver- 
ages (h^^y^^ is equivalent to a Poissonian photon number 



statistics. The change between the two regimes occurs 
when ^{P) changes sign, from negative to positive and 
back, in agreement with the result suggested by the rate 
equation approximation. 

In terms of the phase-space complex argument a, the 
rescaled GS function R{t) changes during the transi- 
tion from being concentrated around the origin into a 
infinitely thin ring-like distribution having a radius of 
yJv{P). For the unsealed GS function i?(s), this corre- 
sponds to a ring distribution whose radius grows in the 
scaling limit faster than its thickness [ss']. This is how 
its value at the origin becomes negligibly small, in accor- 
dance with the discussion in Section IIIII regarding why 
and when is the ansatz of Eq. (fT3|) valid. 



VII. CONCLUSIONS AND DISCUSSION 

The numerical solution of the problem of a two-level 
system in interaction with a cavity mode shows a rather 
sudden change in the behaviour of the steady-state data 
as a function of pumping. Using an analytic approach, we 
reproduce the numerical results and discuss this regime 
change. 

One of the main tools used was the Glauber-Sudarshan 
representation, which was shown to provide the natural 
extension, Eq. (I?7l) . to the recursion relation Eq. ([S]). 
In this perspective it becomes clear that the transition 
observed is related to the vanishing of the GS quasi- 
distribution function R{s) at the origin. 

The laser transition is often represented as the ten- 
dency of the GS function to move away from the origin, 
and this is exactly what takes place in this case too. More 
precisely, it was shown that in the limit k — > 0, 5 — > 
so that = const, the GS function becomes in the 

lasing regime a (5-like distribution concentrated on a ring 
with a pump-dependent radius. In terms of photon num- 
ber statistics this means that the distribution becomes 
Poissonian. 

It is instructive to see the action of the scaling limit 
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directly on the recursion relation Eq. ([8]). Rewritten 
for the rescaled quantities pk = K^Pk its coefficients are 
changed respectively into Ak = A^jn, Bk ~ Bk and 
Ck = KCk- In the scaling limit these coefficients become 
fc-independent 



Ak 



Bk 



P + 7 

-P+7 P+7 



P + 7 
Ck^O, 



in2 



(42) 



with the ratio Ai^/13k converging to !^(P). Then the re- 
cursion relation becomes simply Pk+i ~ v Pk for k \. 
For the values of the pump for which v is negative only 
the trivial solution Pfc = for all A: > 1 is allowed, 
to avoid the appearance of negative expectation values. 
Conversely, when i'{P) is positive a nontrivial solution 
becomes possible. Moreover, according to the above dis- 
cussion, in this case P(0) = and then the additional 
relational = vpo = deduced from Eq. (P7)) also holds. 
Then the solution is pk = i^*^ , in agreement with the state- 
ment that one has i?(t) — > 8[t — v) in the scaling limit 
and the photon statistics becomes purely Poissonian. 

Examining the convergence in Eq. (|42|) it is clear that 
as K becomes smaller and smaller it is the low-fc coeffi- 
cients that get close to their limit first. This is because 
the limit requires the product kn, to be small. There- 
fore the statistics becomes accurately Poissonian for the 
low-index pk first. As k further approaches zero the list 
of nearly-Poissonian expectation values grows longer and 
longer and only in the limit it extends to the whole se- 
quence. In view of these facts g'-^^(O) gets close to unity 
quite early in the process, before the statistics becoming 
entirely Poissonian. As a symptom of a truly coherent 
light, the condition 5^^''(0) = 1 may be, in this sense, 
somewhat premature. 

The scaling limit corresponds to a weak coupling of 
the cavity both to its photon source and to its drain. 
Both coupling parameters involved are present in all laser 
models, so that the limiting procedure can be tested, in 
principle, on any of them, and we have shown its validity 
of the Scully-Lamb model too. As a procedure which is 
successful in two unrelated cases, the scaling limit has a 
good chance of having a wider applicability. 



Appendix A: Differential equation for the 
Glauber-Sudarshan function i?(s) 

The second order differential equation obeyed by -R(s) 
is deduced from Eqs. (|23| as follows. We note first that 



i?„ and T can be eliminated in favour of R and Po the 
first hy Rv — R — Rc and the second using the balance 
condition Eq. (j25|l . As a second step one solves for Rc 
and R'c in favour of R and then identifies the derivative of 
the first with the second. This way Eq. is obtained. 



with the following coefficients 

7 - K + P 

a = K T , 

9 

r 

2P(7 - K + P) 



P 7-K + P / ^ 7-3K + P 



e = 8 + K 



37 - 7k + 3P 
4P(7 - 2k + P) 



7 - K + P 



,7-K + P (7-2k + P)(7-3k + P) 



2g2 



P , ,7-3k + P 
w = 8- - (7 - K + P) 1 



(Al) 



Appendix B: The random injection model 

In this Appendix we apply the scaling limit to the ran- 
dom injection model, in order to show that in this case 
too a sharp transition to Poissonian distribution takes 
place, in the same parameter region as given in the liter- 
ature. 

In the Scully-Lamb laser model [13, HII the cav- 
ity mode is fed by random injection of inverted two-level 
atoms, which exchange energy with the mode through 
Jaynes-Cummings interaction. The time each atom 
spends inside the cavity is a random variable too. The 
losses are simulated by a similar random injection of 
atoms, this time in their lower state. One obtains for 
the time evolution obeyed by the photonic density ma- 
trix 



d 



An 



a/"'"-l + 4„ 



Pn-l.n-l 



A{n + 1) 



A- 1 + f (n + 1)^-" (Bl) 

- Cn Pn.n + C {n+ 1) Pn+l.n+l ■ 

We keep here the notation commonly found in the litera- 
ture [13, m, [131 for the resulting parameters. For B = 
one would obtain the usual competition between the loss 
factor C and the gain (or amplification) factor A. This is 
known to lead to an unstable behaviour for A > C. The 
role of the saturation parameter P > is to limit the 
photon output and to stabilize the solution. A common 
procedure in the analysis of the problem is to expand the 
denominators of Eq. (jBip up to the first order in B, but 
we will not resort to this approximation. 

In terms of the Jaynes-Cummings coupling constant 
one has A and C proportional to and B to g*. Since 
the cavity loss rate C is the same as the parameter k used 
throughout the paper, the limit k — >■ naturally implies 
that g^ ^ too, with their ratio remaining constant. 
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It is easily seen that the steady-state solution of 
Eq. dnil) obeys 



Pn+l,n+l 



Pn,n 



(B2) 



where the notation /i = / BC and (3 — A/B has been 
used. This leads to 



Pn,n 



P" 



(/3 + l)(/3 + 2)...(/3 + n) 



Po,o 



(B3) 



As it was early recognized [17|, the generating function 
of Pn.n is given, up to normalization, by Kummer's con- 
fluent hypergeometric function 



n n 

= Po.Q iFi{l,l3 + l,px) . 



(/3 + l)(/3 + 2)...(/3 + n) 



(B4) 



The value of po,o is determined by the normalization con- 
dition p{x) = 1. The photon number distribution p„^„ is 
connected to the GS function R{s) by Eq. (|22|) which 
implies, with y — 1 ~ x and z ~ y/n 



p{x) 



fc=o 



kl 



■y Pk 



E 

fc=0 



(-1)' k~ 



(B5) 



This shows that p(x) is also the generating function for 
the expectation values pk and of the rescaled ones pk 
when seen as a function of y and z, respectively. There- 
fore p{x) is the appropriate tool for analyzing the scaling 
behaviour of these quantities. According to their defini- 
tion, both /i and f3 scale like 1/g^ so they go to infinity 
as K goes to zero, in such a way that p, ~ up and /3 — njS 
remain finite. 

At this point one can make use of the integral repre- 
sentation of Kummer's function 13 al 



with the notation 



I{z) = I e 

10 



{1 — Kz)up./ K 



{l-uf^'^-^du, (B8) 



and division by /(O) is necessary to ensure the normal- 
ization condition p(0) = p(l) = 1. Rearranging I{z) as 



I{z)^ e-^"''e'^"+'"'"(^-"^'/'"(l-u)-Mw, (B9) 

one brings it into a form familiar from the method of 
steepest descent. The k — asymptotics of /(z) is con- 
trolled by the maximum of ^p{u) — pu -\- P ln(l — u) in the 
interval of integration, since the contribution of the inte- 
grand around the maximum point becomes overwhelm- 
ingly dominant. Depending on the position uq of the 
saddle point, defined by (p'(uo) = 0, i.e. 



Wo 



1-^ 



^ = 1-^ 
p, A 



(BIO) 



one has two regimes, (i) When A < C, uo is negative 
and the maximum of <yj(u) on [0, 1] is found at u = 0. 
(ii) When A> C the saddle point enters the integration 
interval and the maximum of <yj(u) is reached for u — 
1*0 . As a result, and given the normalization condition 
Eq. (jB7l) . one obtains 



p{z) 



ifA^C, 
if ^ > C . 



(Bll) 



In the first case all pk = with the exception of po = 1, 
while in the second one obtains pk — with 



V = puo 



A{A - C) 
BC 



(B12) 



Fi{l,l3 + l,px) ^ 13 e''^"(l - , (B6) 



(B7) 



to express the generation function for pk as 

liz) 



p{z) = p{l - KZ) 



1(0) 



One finds again the same scenario as in Section I VII The 
rescaled expectation values pk obey an exponential k- 
dependence, with a base that changes from zero to a 
positive V at the threshold point. The expression of 
found by the present procedure is the same as given in 
the literature [Ulil, [23|. 
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